Spin- Valley Kondo Effect in Multi-electron Silicon Quantum Dots 
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We study the spin-valley Kondo effect of a silicon quantum dot occupied by TV electrons, with 
TV up to four. We show that the Kondo resonance appears in the M = 1,2,3 Coulomb blockade 
regimes, but not in the TV = 4 one, in contrast to the spin- 1/2 Kondo effect, which only occurs 
at TV = odd. Assuming large orbital level spacings, the energy states of the dot can be simply 
characterized by fourfold spin-valley degrees of freedom. The density of states (DOS) is obtained as 
a function of temperature and applied magnetic field using a fmite-(7 equation-of-motion approach. 
The structure in the DOS can be detected in transport experiments. The Kondo resonance is split 
by the Zeeman splitting and valley splitting for double- and triple-electron Si dots, in a similar 
fashion to single-electron ones. The peak structure and splitting patterns are much richer for the 
spin-valley Kondo effect than for the pure spin Kondo effect. 



I. INTRODUCTION 



Single- and few-electron quantum dots (QDs) coupled 
to leads have been realized in the last several yearsi. 
Their theoretical description is similar to that of a mag- 
netic impurity in a bulk metal, a problem that has been 
studied for decades. It gives rise, among other things, to 
Kondo physics^. The advantage of the QD system is that 
many parameters can be varied by adjusting voltages on 
the electrodes that surround the QD, whereas in the bulk 
these parameters are fixed^. One of these parameters is 
the total number of electrons TV on the QD. Of course 
in many cases some of the electrons on the QD can be 
considered as "core" electrons. This works only when 
the Coulomb interaction or the dot orbital level spacings 
are large and therefore all but one or two of the electrons 
residing on the closest orbital to the Fermi energy on the 
dot can be treated as vacuum. In this single-electron dot 
picture the spin- 1/2 Kondo effect ensues from spin fluc- 
tuation of the dot electron spin coupled to the conduction 
electrons on the leads. As a result, the dot electron spin 
binds to an electron spin on the leads to form a spin sin- 
glet. When temperatures and magnetic field splittings 
are lower than the energy scale characteristic of this spin 
singlet, a narrow zero-bias Kondo resonance crops up in 
the differential conductance as a function of source-drain 
voltage. 

For GaAs QDs^£, the Kondo effect is only observed in 
J\f = odd Coulomb blockade regimes. This is the only 
case in which it is possible for a spin singlet formed by the 
highest-energy electron on the dot and electrons in the 
leads to be the ground state, and then only provided that 
finite orbital spacings render it energetically favorable. 
No Kondo resonance is formed in TV = even Coulomb 
blockade regimes, since the highest-energy dot electrons 
themselves will normally form a spin singlet. Thus a 
characteristic signature of the Kondo effect has always 
been this periodicity of 2 in the variable N '. A second 
characteristic signature is the magnetic field dependence. 
The Kondo resonance, which manifests itself as a peak 



at zero voltage when the conductance is measured as a 
function of source-drain voltage, splits into two peaks 
due to the Zeeman effect. These peaks are separated by 
SV = 2g[i B B. 

For Si QDs, experiments on the Kondo effect have so 
far been able to reach the stage of few-electron dots, but 
not that of single-electron ones. In a few-electron Si doi£, 
a Kondo resonance at zero bias voltage is split into two 
peaks by an applied magnetic field. The value of the 
extracted g factor is about 2.26, slightly larger than the 
expected 2, suggesting a possible contribution from the 
valley degree of freedom. No distinct periodicity has yet 
been observed. 

We note one apparent exception to the periodicity rule 
in GaAs. A recent report^ on the spin-orbital Kondo ef- 
fect in an integer-spin QD has shown a Kondo effect in an 
TV = even Coulomb blockade regimes. Considering two 
orbital levels on the dot, the condition of a spin singlet- 
triplet degeneracy can be artificially achieved by tuning 
the magnetic field and a Kondo resonance emerges at 
a particular magnitude of field. When further increas- 
ing the field, this resonance splits into two peaks at a 
finite bias voltage, and the separation of the two peaks is 
twice the singlet-triplet energy difference. In perpendic- 
ular field, their difference has a much stronger field de- 
pendence than the Zeeman splitting. The smaller-scale 
Zeeman splitting was not observed in the data. The spin 
singlet-triplet Kondo effect involving two orbital levels in 
TV = even regimes resembles in some respects the spin- 
valley Kondo effect in even-electron Si dots, the main 
body of this study, in the sense that the degeneracy does 
not come entirely from spin. 

Silicon is an indirect bandgap semiconductor. The top 
of the valence band lies at the T-point, while the conduc- 
tion band has six degenerate minima along the T-X di- 
rections. The conduction electrons in n-type silicon have 
therefore a sixfold degeneracy corresponding to this val- 
ley degree of freedom. Si QDs are formed in heterostruc- 
tures where the active layer is a thin layer of pure silicon 
sandwiched between layers of Sii^Ge^ alloy, which puts 
the silicon layer in a state of in-plane tensile strain. This 
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raises four of the conduction band minima by a large en- 
ergy (~ 0.1 cV), leaving only a twofold valley degeneracy. 
This in turn is split by small effects that break the mir- 
ror symmetry (reflection through the x-y plane). This 
small (< 1 meV) splitting is enhanced by a perpendicular 
magnetic field. It can be controlled by changing electro- 
static and magnetic confinement. On the other hand, a 
two-dimensional tight-binding model 9 predicts the pos- 
sibility of valley index nonconscrvation during tunneling 
from the leads to the dots, which results in opening up 
additional tunneling channels between an even (odd) val- 
ley state on the leads and an odd (even) valley state on 
the dot, which changes significantly the features of the 
Kondo resonance in single-electron Si QDs. The valley 
mixing is made possible by the rough interfaces that con- 
fine the two dimensional electron gas (2DEG) of strained 
Si, in which the dot-leads system is formed. 

This valley near-degeneracy is a potential source of 
leakage and decoherence in quantum computing schemes 
in which Si QDs serve as the qubits. This is an additional 
reason for trying to understand its consequences. 

In earlier work, we showed that valley degeneracy pro- 
duces a novel Kondo effect in AT = 1 Si QDs^ . We will 
show below that for double-electron Si dots, there is also 
a Kondo resonance, which suffers both the valley and 
Zeeman splittings, in contrast to the Kondo resonance 
affected by the field-dependent singlet-triplet energy dif- 
ference mentioned abovei 

Figure 1 in Ref. characterizes the four spin-valley 
energy levels as a function of magnetic field for a single 
orbital in a Si QD. This energy- level structure draws on 
two experiments^ that observed a valley splitting in Hall 
bars and Quantum Point Contacts (QPCs); in the first 
the valley splitting shows a linear dependence with ap- 
plied field, in the latter the magnitude is measured to 
be about 1 meV. Furthermore, a finite zero-field valley 
splitting was observed in both experiments, in which at 
zero field the splitting in Hall bars is about 1.5±0.6 ^teV, 
much smaller than in QPCs . The difference is ascribed 
to interface disorde r 10 ! 11 . This level structure will again 
be used later to define the dot energy levels in Sec. IIVI 

Our aim in this paper is to describe the Kondo reso- 
nance^) in Si QDs for TV > 1. The Kondo effect itself 
results from a ground state in a which the dot spin and 
valley states are mixed with the lead states to form a 
singlet ground state. The resonance in transport occurs 
because the resulting wave function has weight at or near 
the Fermi energies of both leads, leading to a zero-bias 
or near-zero-bias anomaly. The anomaly can be shifted 
and split by a magnetic field. For AT = 1, we showed 
that the zero-bias resonance is split by both the valley 
and Zeeman splittings so that there are more split peaks 
in non- linear I — V characteristics than in the spin-1/2 
Kondo effect^. To complete the whole picture, we con- 
tinue to investigate whether there is a spin-valley Kondo 
effect in the AT = 2, 3, 4 Coulomb blockade regimes and if 
so, how the splittings occur in each regime. We use again 
an equation-of-motion approach to obtain the interacting 
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FIG. 1: When the Coulomb interaction is large, the electron 
occupation number AT is a good quantum number at weak 
couplings. Shown here are the six and four configurations of 
low-energy spin-valley dot states for the (a) TV = 2 and (b) 
TV = 3 regimes, respectively. The spin-valley Kondo effect 
is related to the TV-conserved many-body transitions where 
these dot states interact with the lead states. Note that the 
relative size of the valley and Zeeman splittings in the cartoon 
figures does not necessarily reflect every tested sample since 
the valley splitting is sample-dependent. 



DOS that will be shown below to directly reflect the dif- 
ferential conductance. We assume orbital level spacings 
to be larger than a few times the dot Coulomb interac- 
tion U, so that the orbitals are well separated in energy. 
As a result, it is sufficient to only consider the single or- 
bital closest to the Fermi energy. This condition could 
be satisfied in a small dot size of ~10 nm. Then, since 
the Coulomb interaction is large, the charge fluctuation 
is small enough for the electron occupation number TV to 
be a well-defined quantum number between the Coulomb 
blockade peaks. Then what needs be taken into account 
is the low-energy electron configurations for each TV. Fig- 
ure[T]demonstrates the configurations of these spin- valley 
energy states for TV = 2, 3 that interact with the lead 
states. 

In the next section we describe the model and our cal- 
culation method. The results for the Kondo temperature 
follow in the section after that, and then we plot and dis- 
cuss the DOS for various cases. Finally we summarize 
and give some conclusions. 
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II. EQUATION-OF-MOTION APPROACH 

The equation-of-motion (EOM) approach has proven 
to be a good tool for investigating the density of states 
(DOS) in the Kondo effect^. Since this quantity is the 
one that interests here, we shall employ this method. In 
particular it is able to handle arbitrary orbital structure 
and particle number, as well as finite temperature and 
the presence of a magnetic field. The basic technical 



details are given in Ref. 0, so here we describe only those 
additional features that are required to treat the J\f > 1 
case. 

A Hamiltonian that describes a system consisting of 
the single-particle energy levels of the leads, the dot, the 
tunneling matrix couplings that connect the levels of the 
leads and the dot, as well as the Coulomb interaction 
between electrons on the dot is the Anderson impurity 
model, expressed by 
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Here the spin-1/2 index a <G {f, J.}, and the valley in- 
dex m £ {e(even), o(odd)}. Even and odd denoted the 
two valley states, fh is the opposite of the valley index 
to. The operator cl kma (cih ma ) creates (annihilates) an 
electron with an energy Sk in the i lead, i £ L,R, while 
the operator f ma (f mC r) creates (annihilates) an electron 
with an energy e ma on the QD, connected to the leads 
by Hamiltonian intravalley coupling Vo,ik and intervalley 
coupling Vx,ik- For perfect interfaces at the boundaries 
of the well, it could happen that Vx,ik = 0, i.e., that the 
valley index is conserved in tunneling. For real dots, we 
have shown that Vo.ik and Vx,ik arc likely to be the same 
order of magnitude^. We assume that Vo(x)ik does not 
depend on the spin index a. U is the Coulomb interac- 
tion on the dot and is assumed to be independent of the 
valley index. 

The Kondo effect can be observed by measuring the 
current / and likewise the differential conductance G 
as a function of source-drain voltage V s d- Theoretically 
the differential conductance G = dI/dV s d is given by 
differentiating the generalized Landauer formula, given 
in Ref. ||J The differential conductance is approxi- 
mately proportional to the interacting DOS — DOS 
= — lm[Q mcr (eV s d)]/'!T — given the assumption of an ini- 
tially flat noninteracting DOS in the leads. Although the 
lead DOS and tunneling matrix elements vary with ap- 
plied voltages, it is usually true that the variations are 
slow compared with the sharp Kondo resonance struc- 
tures. When this is true, understanding the dot DOS is 
sufficient to identify the fine structure in the conductance 
near zero bias. 

Thus we need to compute Qma{w), the retarded 
Green's function: 

G ma {w) = {{fma, f La)) 

= -i / e iat {{f m(7 {t)Jl a {0)})dt. 
Jo 



where a = w + iS. We compute the equation of motion 
for O m a{w) in the frequency domain: 

w((A,B)) = ({A,B}) + (({A,H],B)) 
= ({A,B}) + ((A,[H,B])). 

By applying the above equation of motion to a Green's 
function, we obtain higher-order Green's functions on the 
right-hand side of the equation, which we further expand 
by repeating the same procedure until all second-order 
(0(V 2 )) contributions are preserved after the decoupling 
scheme^. After some tedious but straightforward calcu- 
lations, we acquire equations of motion that couple the 
Green's functions <7 m<T (to), {{fma, fLa))> ii n ifma, f mtT )), 
{{ni'fmaj ma )), {{nWjfmcf^}}, {(ni>nj>f ma , p m<r )), 

{{nifljUpfmaJ^}}, {{ni'nj'Tlp'ffncfhc,}) ({IJiP} + 

ma and {/', j', p'} ^ ma are shorthand notations for both 
to and a indices) that describe single, double, triple and 
quadruple occupancies, respectively. These Green's func- 
tions suffice to describe fourfold degenerate Si QDs that 
can host up to four electrons. They are nothing but lin- 
early coupled matrix elements spanned by spin and valley 
quantum numbers, and they can be solved for by lin- 
ear diagonalization in terms of their coefficients: second- 
order perturbation terms, integral functions, occupa- 
tion numbers (n mcr ), (nin mrT ), (ninjn mtT ), and expec- 
tation values (f m(T fma), {ni" f ma fma), ("("^"/mff/mcr) 

ma i ma) . It is noteworthy that some perturba- 
tion terms and the integral functions are logarithmically 
divergent at the Fermi energy, thereby giving rise to a 
zero-bias anomaly, and have to be treated carefully in 
the Kondo regime T < Tk {Tk is the Kondo temper- 
ature that we will define later). Details are parallel to 
those in Ref. so that we do not give them here, except 
to note the following: 

First, we assume a flat and symmetric noninteract- 
ing DOS in the source and drain, so Vo(x),ik = ^b(x)- 
Since the valley index is not conserved, it is convenient 
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Vcos<p, Vx = Vsm< 



We define 



V§ + V x and 2V V X = PV 2 with f3 = sin 20 and 



to introduce Vo 

v 2 = v 2 + v x 

< f3 < 1. Parameter /? gauges the extent of valley 
index nonconservation. In other words, /3 = implies 
valley index conservation, reflecting perfectly smooth in- 
terfaces of the 2DEG, whereas (3 = 1 the maximal valley 
mixing. 

To compute numerical results, we use the following in- 
tegral 



dw — — = -W(-± 

_ D w — W 1 ±10 



w — w 



2iriT 



D + w in 



where the parameter D is the conduction half-bandwidth, 
and /fd the Fermi function, ^(z) is the digamma func- 
tion that asymptotically behaves as Inz as \z\ 3> 1. This 
logarithmic divergence produces a low energy scale in the 
Kondo regime, the domain of our interest. This scale is 
defined as the Kondo temperature. We also assume the 
self-energy term 

V 2 

s oM = V — r ^ -«r 

W - £ k + 10 

where V = irV 2 /D. This approximation is valid near the 
Fermi energy which is the region of most experimental 
interest. 

The integral functions come from the correlation func- 
tions of the dot and leads after decoupling the higher- 
order Green's functions. By simple transformation these 
correlation functions can be rewritten in terms of inte- 
gral functions over the Green's functions shown above, 
and the set of equations of motion terminates after the 
decoupling^. To compute the integral functions whose 
integrands contain the Green's functions, first we use the 
approximation adopted by Lacroi»i^ and V. Kashcheyevs 
et alr^ who assume that the integral functions are dom- 
inated by the singularity and only strongly affect the 
region around the Fermi energy. As a result , we can 
approximate, for instance, the integral function 



dw 



,f F D{w')g; n(T {w') 



TQ* ma {e F ) / dw 



JfdW) 



w — w' 



(1) 

Likewise for other integral functions. 

This leads to a set of coupled integral equations for 
the exact Green's functions since they are functions of 
integrals of themselves. Our strategy here is to make a 
guess of the expected structures of the Green's functions, 
substitute them for the exact ones in all the integral func- 
tions, and iterate to self-consistency. To illustrate, we re- 
place the Green's function Q ma (w) in Eq. [T]by Gmaiw) 
in Reffl5l as the expected Green's function. If we assume 
T/U <C 1, we find Q ma (w) consistent with the features 



of Gma(w), and equivalent to GmV{w) at high temper- 
atures where the integral functions are small and can 
be disregarded. Moreover, to mimic the Green's func- 
tions that feature a broad peak of width ~ T centered 
around the discrete bare dot energy levels, each propaga- 
tor [w — s mrT — dU]^ 1 in Gma(w) is given a finite spectral 



(at). 



width Td where d = 0, 1, 2, 3. T^s take into account only 
the self-energy terms Eo(w)'s in each propagator. Thus 
we assign 

r = r, r 1 = 3r, r 2 = 5r, r 3 = 7r. 

This increase of spectral widths mimics real physical sys- 
tems where the potential barrier gradually opens up at 
high energy. Despite the fact that T^'s are not real spec- 
tral widths, the errors contribute only 0(T 2 ) to the inte- 
gral functions that already have a prefactor T(c.f. Eq.[T]). 

We shall operate under the assumption that T/U <C 1, 
yet we are interested in the Kondo regime. Hence this ap- 
proximation for the integral functions might not be valid. 
However, as argued by CzycholU^, the approximated in- 
tegral functions should not deviate very much from the 
true ones, so long as the temperature is not too far be- 
low the Kondo temperature and higher than second-order 
contributions can be legitimately neglected. 

On the other hand the occupation numbers and other 
expectation values are computed by integration over the 
Green's functions. For instance 



(nin ma ) = 

ifmafmij) = 
i n l" fmafmc) = 



f 



IT 



7i 



dwf FD (w)lmQ ma (w) . 
dwf FD (w)lm((nif ma , fl lrJ )). 

i J dwf FD (w)lm((frn a , f} na )). 

^ J dwf FD (w)Im((ni"f ma , f} na ))- 



Now that all the perturbation terms, integral functions, 
the occupation numbers, and other expectation values in 
the Green's functions have been accounted for, we pro- 
ceed to iterate to self-consistency. 

In Sec. IIII1 we show the Kondo temperatures for JV = 
1,2,3 regimes as a function of finite Coulomb interaction 
U, dot energy level s ma , coupling strength T, and param- 
eter /3. We also present numerical results of the DOS as a 
function of energy, in which a Kondo resonance is found 
at zero bias voltage and splits when applying a magnetic 
field, as demonstrated in Sec. HVl 



III. 



THE (^-DEPENDENT KONDO 
TEMPERATURES 



It is enlightening to start first with a simple case: de- 
generate spin-valley states. Although it is unlikely to 
observe a degenerate spin-valley Kondo effect in Si dots 
since the valley splitting is not zero even in the absence 
of external field, it is still useful for us to first derive the 
Kondo temperatures from the solutions near the Fermi 
energy of the real parts of the denominators in Q ma (w)—. 
Assuming D > |e mo -|,f7 3> Tr-'s and at zero tempera- 
ture, we obtain 
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T± (13) ~ n* oxo 7r ( £ "^- £ ^)( £ ^ + ?7 - £ ^) 

K2 ^j ~ exp (3 ± ^p^^g^ _ £p) + (£mCT + 2U _ £p)] 

T ± 7r(£ mg . + E7 - e f )(e mg + 2U- e F ){e ma + 3U - e F ) 

KaW ~ cxp (3 ± ^ )rc/[(£m<7 + ^ _ £f) + 4/3(£mCT + 3t/ _ £f)] 



Where D* = D\2e mrT + U - e F \/\D + 2e mtJ + U\. 
T^(P), T| 2 (/3), T± 3 ((3) are the /3-dependent Kondo 
temperatures for single, double, and triple occupancies, 
respectively. We have already shown 9 - that there are two 
Kondo temperatures — equivalent to T^ 1 (/3) in the infi- 
nite U limit — in single-electron Si dots if the valley in- 
dex is not conserved (j3 > 0), and the parameter /? greatly 
enhances one Kondo temperature while suppressing the 
other in a similar fashion. As a result, T^ 1 (/3), the largest 
of the two, dominates the screening of the dots. Here we 
show the same with double-electron and triple-electron Si 
dots. But when e ma + 3U < e F , the spin- valley Kondo 
effect disappears since all energy levels are occupied and 
inelastic transitions are prohibited. A keen observation 
will give that these T±(/3)'s are not unrelated. Indeed, 
they obey a simple relation 



1 



lnT± 2 (/?)/£* 



\nT^(P)/D* 



(2) 



In other words, the Kondo temperatures in the Af regime 
are influenced by the ones in the neighboring Af — 1 and 
Af + 1 regimes. A similar logarithmic relation was found 
in Rcf . [l6|, where although a Kondo effect is produced in 
Af = even regimes, it belongs to a two-level QD. Conse- 
quently it is gcncrically different from a fourfold degen- 
erate Si QD. 

The larger Kondo temperatures Tj£(/?)'s express the 
energy scales where the spin- valley Kondo effect can be 
observed in the Kondo regime, namely, at temperatures 
lower than T^(/3)'s. The EOM approach produces a 
finite-?/ Kondo temperature T J ^ 1 (/3 = 0) that is different 
by a prefactor 3/4 in the exponent from the one using a 
scaling theory^. The discrepancy of this approach was 
claimed to be due to neglect of higher-order corrections^. 
It should therefore be stressed that Tj£(/3)'s may only 
show a qualitative dependence on the model parameters. 
Despite this defect, T^ 1 (/3)'s computed from the EOM 
approach scale logarithmically the zero-bias resonance in 
single-electron Si QDs£, thereby portraying a Kondo-likc 
effect, as will be expected the same with T± 2 (/3) and 
T| 3 (J3). The larger T+ (/?) 's will serve as energy scales for 
the demonstration of the spin- valley Kondo resonances in 
the next section. 

It is noteworthy that at zero temperature the above 
solutions exist except when e ma = (ep — U)/2. This par- 
ticular condition leads to a particle-hole symmetry. The 



EOM approach fails to produce at this point an energy 
scale that governs the low temperature behavior. This 
scenario requires a special mathematical treatmen t 14 ' 19 . 
Regardless of its intrinsic interest, this case is not our 
main concern in this paper. 



IV. DENSITY OF STATES 

In finite magnetic field, the dot energy levels arc given 
by e ma =e d + (A/2 + ^ v B)(S m ^ - S m . ) + gfjL B B(8 a ^ - 
S ay i), where Ed is the dot bare energy level and B is the 
applied magnetic field. A is the zero-field valley splitting 
and n v is the constant valley splitting slope. Note that A 
and fj, v constants are sample-dependent and may vary in 
actual experiments. In the following, we have taken fj, v = 
0.1 mcV/T, slightly smaller than g^ B = 0.114 mcV/T 
where g = 2. We have also chosen the parameters T = 0.2 
meV, D = 200r, and U = 40I\ which, though bringing 
forth small Kondo temperatures, are favorable for the 
purpose of illustration. 



Degenerate spin-valley Kondo effect with valley 
index conservation 



We have obtained in the previous section the /in- 
dependent Kondo temperatures. It is now instructive 
to apply them numerically as a temperature scale to the 
Kondo resonance. Our task is to plot the DOS, each 
time tuning the bare energy level Ed in order to shift 
the Fermi energy ep (which we set to be zero) inside 
the Af = 1,2,3 Coulomb blockade regimes, mimicking 
experimental manipulation of the gate voltage over the 
dot energy levels. Let us for a moment consider valley 
index conservation and full spin- valley degeneracy by dis- 
regarding the expected zero-field valley splitting. Figure 
[5] displays three plots of the DOS near the Fermi energy 
in the Af — 1,2,3 regimes. Not surprisingly, for odd 
Af (1,3), a narrow Kondo resonance at zero bias voltage 
can be seen in Fig. 2(a)|and Fig 2(c) Another side peak 
around 2 meV in Fig. |2(a) comes from the process that 
two electrons are depleted from the dot simultaneously 
and therefore is energetically disfavored at large U. It 
is interesting to note that the positions and shapes of 
the Kondo resonances for Af = 1,3 reflect a particle-hole 
symmetry. Also noteworthy is that in Af = odd regimes 
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2: The DOS's with approximately 2(a) < 
2(c) three electrons on the dot. Here the valley and Zeeman 
splittings are assumed to be zero. The valley index is consid- 
ered conserved. In these three regimes there is a Kondo reso- 
nance at the Fermi energy, shown in the insets. To plot them 
we use |2(a)| e d = -15r, T/T K i = 0.2 with T K \ ~ lxl" 
e d = -55F, T/T K2 = 0.2 with T K2 ~3x 10" 4 
e d = -105r, T/Tks = 0.4 with T K 2 ~ 1.6 X 10~ 4 meV. 
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The self-consistently computed occupation numbers are 2(a) 
< n ma >= 0.2505 « 1/4; [2(b)] < rw >= 0.4831 « lj2 
[2(c7| < n mrT >= 0.7312 



the dot displays itself as a spin-1/2 magnetic impurity 
and is fully Kondo-screened by the conduction electrons 
on the leads. The pseudospin possessed by the valley 
degree of freedom behaves similarly. 

On the other hand, a more symmetrical zero-bias reso- 
nance in the W = 2 regime than the two in the N = 1, 3 
regimes is found in Fig. 2(b) This unexpected resonance 
is attributed to the valley degree of freedom that provides 
valley states for additional inelastic transitions to occur: 
for example, spin flips can occur even when the ground 
state is a singlet, if the flip is accompanied by a change in 
valley index. In fact, single-electron and double-electron 
Si QDs belong to the types of fully-screened, and under- 
screened Kondo effect, respectively. 

A somewhat similar phenomenon of the Kondo ef- 
fect can occur in a two-level quantum dot with orbital 
degeneracy^. By introducing an interpolative pertur- 
bative approach the Kondo effect is obtained in the 
J\f = 1,2,3 regimes. However, this approach does not 
give an accurate estimate of the Kondo temperatures as 
a function of model parameters and only provides a qual- 
itative thermal behavior. The EOM approach we have 
chosen here acquires similar structures for the DOS and 
seems more favorable to render better qualitative (if not 
quantitative) Kondo temperatures and Kondo resonance 
structures. 



B. Field dependence and effect of valley index 
nonconservation for the J\f = 2 regime 

The splitting of the Kondo resonance is instructive be- 
cause it disentangles clearly the interplay of all the par- 
ticipating low energy states and their magnetic field de- 
pendences. The peak splittings are due to the lifting of 
valley and spin degeneracies. Here we demonstrate the 
field-dependent peak splittings while taking into consid- 
eration the effect of the valley index in two extreme cases: 
conservation (/3 = 0) and full nonconservation ((3 = 1) of 
valley index. We shall particularly concentrate on the 
unusual Kondo peak in the J\f = 2 regime. 

Consider first (3 — 0. Fig. [3] demonstrates the peak 
splitting structure of double-electron Si QDs. The zero- 
bias peak splits into three peaks in red curve when the 
zero-field valley splitting A ^ and the magnetic field 
B = 0. Among these three peaks, the central splits into 
two peaks by the Zeeman splitting when B ^ 0, since 
it originates from a spin Kondo effect, whereas each side 
peak splits further into three; see the blue curve. Each 
peak correspond to a many-body transition. The argu- 
ments parallel those for M = 1—: the spin-valley Kondo 
effect comes from spin-flip, intcrvalley, and intravallcy 
inelastic transitions among the dot states in Fig. 1(a) 
interacting with the lead states. The valley and Zee- 
man splittings break the degeneracy of these dot states 
to produce additional peaks. We show the transition 
corresponding to each peak in the schematic diagrams of 



Fig. 3(b) 
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FIG. 3: (Color online) The DOS for Af = 2 is plotted 
with P = (valley index conservation), the zero-field val- 
ley splitting A = 12 Tki and a magnetic field B = (red 
3 x 1CT 3 T (blue curve). The Kondo temperature 



curve 

T K i ~ 3 x 10 
for Fig. [2(b) 



4 meV. T/T K 2 = 0.2, e d = -55F, the same 
The split peaks and their field dependences 
parallel those of single-electron Si QDs in Ref. 0, with similar 
many-body transitions in the schematic diagrams a, b, c, d 
producing peaks a, b, c, d, respectively. 



FIG. 4: (Color online) The DOS for JV" = 2 is plotted with 
(3 = 1, the zero-field valley splitting A = 10 T+ 2 (/3 = 1) 
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T (green 
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and a magnetic field B = (red curve), 5 x 
curve). The Kondo temperature T^ 2 (f3 = 1) ~ 5 X 10~ a meV 
T/T+ 2 (/3 = 1) = 0.2, e d = -55r. Note that there emerges 
an extra zero-bias peak due to the effect of valley index non- 
conservation, compared with the peak splitting structure in 
Fig(3] where the valley index is conserved. It introduces an- 
other transition that involves an electron hopping between 
opposite valleys, as shown in the schematic diagram. 
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cases are rather similar. This is due to particle-hole 
symmetry: instead of electrons, it is the holes that tun- 
nel in and out of the dot. The four spin- valley energy 
states interac ting w ith the lead states for Af = 3 case are 
shown in Fig. 1(b) With these four states gradually sep- 



arated by the valley and Zccman splittings, the inelastic 
co-tunnelings produce a similar peak splitting structure 
to that for A^ = 1, a case which has been exhaustively 
treated in Ref. To avoid redundancy, we omit the plot 
of its DOS. 



CONCLUSIONS 



At j3 — 1, we have maximal nonconservation of the 
valley index. For this case, Figured] shows an additional 
zero-bias peak that represents a bound state of opposite 
valleys, which is therefore a a manifestation of the pure 
valley Kondo effect, apart from the other eight peaks 
already identified in Fig. [3] The transitions that gen- 
erate this peak take the following course: a conduction 
electron at the odd (even) dot state tunnels out to the 
odd (even) lead state through the intravalley coupling 
Vo, while another electron in the even (odd) lead state 
tunnels into the odd (even) dot state through the inter- 
valley coupling Vx ; see the schematic in Fig. [4] This 
peak height increases along with the intcrvallcy coupling 
Vx j or /?, which would therefore provide an experimental 
signature of valley index nonconservation. 

As already seen in Fig. [21 the N = 1 and TV = 3 



We summarize our results as follows: (a) The spin- 
valley Kondo effect appears as expected in the Af = 1 
and Af = 3 regimes, and unexpectedly also in the Af = 2 
regime, but not in the Af = 4 regime, thus yielding in 
general a Kondo effect unless Af is divisible by 4. This 
contrasts with the spin Kondo effect which appears only 
at odd Af. (b) Figure [2] shows an asymmetrical struc- 
ture astride the Fermi energy in the position and shape 
of the zero-bias Kondo peaks in the Af = 1 and Af = 3 
regimes. In the Af = 2 regime the peak shape is more 
symmetrical, (c) By applying a magnetic field, the peak 
splittings in the Af = 2 regimes resemble that in the 
Af = 1 regime. We are able to to attribute each peak 
to its corresponding inelastic many-body transition (co- 
tunneling); see the cartoon schematics in Fig. [3l The 
rich level structure gives rise to a rich pattern of peaks. 

We expect that these many-body signatures of the val- 
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ley degree of freedom in Si will be observed in future 
experiments. 
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